Fermi liquid behavior of the normal phase of a strongly interacting gas of cold atoms 



O 

(N 

00 



C/2 



o 

> 



X 



S. Nascimbene^*, N. Navon\ S. Pilati^^ F. Chevy\ S. Giorgini^, A. Georges^'^, C. Salomon^ 
^ Laboratoire Kastler Brossel, CNRS, UPMC, Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris, France 
^ Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 
^ Dipartimento di Fisica, Universita di Trento and INO-CNR BEC Center, 1-38050 Povo, Trento, Italy 
* Centre de Physique Theorique, CNRS, Ecole Poly technique, route de Saclay, 91128 Palaiseau Cedex, France 
^ College de France, 11 place Marcellin Berthelot, 75005 Paris, France 
(Dated: June 9, 2011) 

We measure the magnetic susceptibility of a Fermi gas with tunable interactions in the low 
temperature limit and compare it to quantum Monte Carlo calculations. Experiment and theory 
are in excellent agreement and fully compatible with the Landau theory of Fermi liquids. We show 
that these measurements shed new light on the nature of the excitations of the normal phase of a 
strongly interacting Fermi gas. 
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In 1956 Landau developed an elegant description of 
interacting Fermi systems at low temperature relying 
on the existence of long-lived quasiparticles. While this 
Fermi liquid theory (FLT) describes well Helium 3 and 
many solid-state materials above the superfluid temper- 
ature, there exist notable exceptions such as underdoped 
cuprates [l|, where despite tremendous theoretical and 
experimental efforts, the nature of the normal phase is 
not yet understood. Similarly to high critical tempera- 
ture superconductors, the properties of the normal phase 
of strongly correlated atomic fermionic gases and the na- 
ture of its excitations are still debated. This issue was 
addressed recently for spin-balanced gases above the su- 
perfluid transition, through the measurement of equa- 
tions of state 0, y-IlL study of the single-particle 
excitation spectrum or of spin fluctuations [lo| . 

On the one hand, recent photoemission spectroscopy ex- 
periments near the critical temperature were interpreted 
using a pseudogap model On the other hand, mea- 
surement of the temperature dependence of the specific 
heat displayed a linear behavior compatible with Fermi 
liquid's prediction 2]. All these experimental probes give 
access to the properties of the normal phase of the un- 
polarized normal phase above the critical temperature 
Tc- This limitation can be overcome by stabilizing the 
normal state at T < Tc by i mposing a spin population 
imbalance in the trapped gas llHl3[ and extrapolating 
its properties to zero imbalance. Previous works focused 
on the highly-polarized limit where minority atoms be- 
have as impurities: n2 ^ jt-i, where rii is the density for 
species i 0, 0, Ei " 21 j . Here, we interpret the spin imbal- 
ance as the application of an effective magnetic field to 
the unpolarized normal gas at very low temperature and 
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using a combination of Monte Carlo simulations and ex- 
perimental results, we extract from the equation of state 
the magnetic spin response of the normal phase in the 
limit T <^Tc. We show that our results are compatible 
with a Fermi Liquid description of the normal phase, and 
we extract the Fermi liquid parameters in the universal 
unitary limit where scattering length is infinite. The re- 
lationship between these parameters and the properties 
of low lying excitations of the system allow us to quanti- 
tatively interpret spectroscopic data from Q . 

The polarization dependence of the energy E of 
the system directly reflects the presence of spin-singlet 
dimers in the sample. Indeed, the presence of a gap in the 
spin excitation spectrum implies a linear dependence of 
the energy E with polarization p — [Ni — N2)/{Ni -f iV2) 
at low temperature, and hence a zero spin susceptibility. 
We have performed quantum Monte Carlo simulations of 
the partially polarized Fermi gas at T = in the BEC- 
BCS crossover. We make use of the fixed-node diffusion 
Monte Carlo method that was employed in earlier stud- 
ies of polarized Fermi gases 
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The state of the 
system is forced to be in the normal phase by imposing 
the nodal surface of a many-body wave function incom- 
patible with off-diagonal long-range order. A simple way 
to implement this requirement is by choosing the trial 
function of the Jastrow-Slater form 



(1) 



where R = (ri, r^) is the spatial configuration vector 
of the N particles and D denotes the Slater determinant 
of plane waves in a cubic box of size L with periodic 
boundary conditions. The positive Jastrow correlation 
term f(r) is determined as described in Ref. [l5| : at short 
distances it corresponds to the lowest-energy solution of 
the two-body problem, while it satisfies the boundary 
condition on its derivative fir = L/2) — 0. 

The results for the canonical equation of state 
E{Ni,N2) are shown in Fig[TJ They are well fitted by 
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FIG. 1: Canonical equation of state of a two-component Fermi 
gas calculated using quantum Monte Carlo simulation, for 
l/kpa = -1.5, -1, -0.6, -0.2, 0, 0.2, 0.4, 0.5 (from top to bot- 
tom). The solid lines are fits of the low-polarization data with 
equation ((2)1. Inset: Extracted values of the susceptibility x 
as a function of l/kpa. The dashed red line is the result of a 
perturbation expansion valid up to order (kpaf' ■ 
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holding for a spin polarizable system at low temperature, 
where both ^jv and the dimensionless spin susceptibility 
X (in units of the susceptibility of an ideal Fermi gas 
■in/2EF) depend on l/kpa, where kp = {iir'^nf/^. The 
Monte Carlo method indicates the absence of spin gap, 
and thus of preformed molecules in the normal phase 
for l/kpa < 0.5. Note that the extracted values of x 
reported in the inset of FiglDshow a rapid drop for pos- 
itive values of a when entering the BEC side of the Fes- 
hbach resonance. A likely explanation is the binding of 
fermions into spin-singlet pairs for some positive value of 
the interaction strength l/kpa. Monte Carlo calculations 
for values of l/kpa > 0.7 show that E{p) is indeed lin- 
ear rather than quadratic in p, indicating the emergence 
of a gap. However, pairing fluctuations play a major 
role for such values of the coupling and the nodal surface 
of the Jastrow-Slater state (see Auxiliary Materials) is 
no longer sufficient to enforce the normal phase. This 
behavior is reminiscent of the pairing transition investi- 
gated in the framework of BCS theory [2^, as well as in 
the normal phase of the attractive Hubbard model, ex- 
trapolated to a temperature range below the superfluid 
transition 23, 2^ , while in our work the extrapolation is 
made towards a small spin imbalance. 

We now compare these simulations with the grand- 
canonical equation of state (EoS) of a homogeneous sys- 
tem obtained experimentally in ref. 0, 0] ■ We prepare 
a deeply degenerate mixture of the two lowest internal 
states of ^Li, held in a cylindrically symmetric hybrid 
optical/magnetic trap, of radial (axial) frequency (w^, 
respectively). The bias magnetic field Bq is chosen be- 



FIG. 2: Thermodynamic function h{b) measured at different 
magnetic fields Bo = 871, 834, 822 G. The blue lines corre- 
spond to the superfluid equation of state /is ((5) measured in 
[j. The red line is a linear fit of the data in the normal phase, 
b > be. The dashed line indicates the superfluid/normal phase 
transition (6 = 6c). 



tween 822 G and 981 G, allowing to tune the strength of 
interactions —1 < l/kpa < 0.2. The final atom number 
is 2 to 10 X 10* atoms per spin state, and the gas tem- 
perature is smaller than 0.06 Tp, as measured from the 
fully-polarized wings of a trapped gas [Hj. From dimen- 
sional analysis, the EoS of a spin-imbalanced Fermi gas 
can be written as 



P(^i,^2,a) =^ Poi^J,)h[S 



y/2miia 



where /i = (/ii + /i2)/2 is the mean chemical potential 
and Po{fJ-) is the pressure of a non-interacting unpolar- 
ized Fermi gas. i5 is a grand-canonical analog of the inter- 
action parameter l/kpa, and 6 is a dimensionless number 
proportional to the 'spin-polarizing field' /Zi — IJ.2- 

At all values of the scattering length addressed in this 
work, the equation of state exhibits a clear discontinuity 
of its derivative at the critical field bc{5) (See Fig. [5]), 
indicating a first-order phase transition from a superfluid 
state for 6 < 6c to a normal state for b > be where h 
is linear in 6^. 0, [l3|- The equation of state of the 
superfluid phase has been discussed in a previous work 
Q and we focus here on the properties of the normal 
phase. We write: 



h{S,b) = hj,{S)[l + ^x''''{S)b' 



0{b^) 
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hf^{5) is the grand-canonical equation of state in the nor- 
mal state, extrapolated to a spin-symmetric configura- 
tion. x^^((^) is a grand-canonical magnetic susceptibil- 
ity. For an ideal two-component Fermi gas, the functions 
and x'^'-^ are equal to 1. Fitting our data in the nor- 
mal phase with ([3]), we obtain the parameters /iAr((5) and 
X*^^(5) in the BEC-BCS crossover shown in Figl3] where 
we compare their values to the predictions of the Monte 
Carlo simulations. To this end, we fit the dependence 
with l/kpa of the parameters ^at and x determined by 
Monte Carlo simulations, and perform a Legendre trans- 
form to obtain the grand-canonical EoS hj^i [5) of the nor- 
mal phase and magnetic susceptibility x'^'^i^) measured 
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FIG. 3: Fermi liquid equation of state extrapolated to a spin- 
symmetric configuration /ijv((S). The black dots are the exper- 
imental data, and the red line is calculated from the Monte 
Carlo data. Inset: Grand-canonical susceptibility x'^'^i^) of 
a Fermi gas in the BEC-BCS crossover. 



experimentally. In the investigated parameter range, the 
agreement between theory and experiment is excellent. 
We also remark that our value for the susceptibility of 
the normal phase at unitarity is about twice larger than 
the value measured in on a gas with a 35% conden- 
sate fraction, confirming a significant suppression of the 
spin susceptibility in the superfluid phase . 

Our findings demonstrate that for 1/kpa < 0.5, the 
spin excitations of the system are not gapped in the 
normal phase which therefore does not support "true" 
molecules. However, a certain class of theories predicts 
a reminiscence of this gap in the form of a dip in the 
density of states over a range A* around the Fermi level 
[26| . A* is often called the pseudogap, in relationship to 
some features of high-critical temperature superconduc- 
tors. These theories predict a departure of E{p) from its 
quadratic behavior when the Fermi levels of the two spin 
species reach the edges of the dip, /i2 — /^i — A*, (see 
Auxiliary Materials). The absence of such an anomaly 
in Fig. I1I2[ and in the whole range —1 < 1/kpa < 0.5 
thus suggests that the dip is either extremely narrow or 
very broad: the density of state remains flat over the 
range of polarizations and interaction strength studied 
in our work. For instance, at unitarity this range covers 

< &^ < 3. If a sizeable dip existed, then its width can- 
not be smaller than ~ (/ii + /i2)\/3 — lAEp where we 
have used the unitary equation of state, /z = 0.41 Ep 

01 • Such a large pseudogap is not compatible with 
the photoemission data of [9|(See below). Furthermore, 
we would expect on physical grounds that A* becomes 
smaller on the BCS side of the resonance. This is ob- 
served neither in the experimental data of Figl2] nor in 
the Quantum Monte Carlo results of Fig[T] 

On the contrary. Landau's theory of Fermi liquids is 
fully compatible with our observations. This theory as- 
sumes the existence of long-lived fermionic excitations 
above the Fermi surface. Combining the measurement of 



the low-temperature compressibility k and specific heat 
Cy of with the data presented here, we can fully 
characterize the parameters of the theory at the uni- 
tary limit. From the magnetic response of the T = 
gas, we obtain here its magnetic susceptibility and an- 
other determination of k. The two determinations of 
K coincide within 5%, showing that the two approaches 
indeed probe the same Fermi liquid. From this set of 
thermodynamic quantities we derive, according to Lan- 
dau's Fermi liquid theory, a complete characterization 
of the low-lying excitations of the unitary gas: besides 
their effective mass m* = 1.13 rri and Landau param- 
eters = -0.42, = 0.39 found in 0, we re- 
—0.40 and obtain the new parameter 
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Fg- = m*/mx{0)-^ - 1 = 1.1(1). Note that Eg- > 
corresponds to magnetic correlations which do favor the 
singlet configuration. 

We can finally test FLT on the single-particle photoe- 
mission spectrum obtained at the unitary limit and at 
the onset of superfluidity from ref. . The experimen- 
tal signal A{k,Lu) is directly proportional to the spectral 
function A{k, uj — ^) averaged over the trap that we es- 
timate using the following procedure: In the vicinity of 
the Fermi surface, the dispersion relation of the Fermi 
liquid quasi-particles reads fiuik = Ai + ?i^(fc^ — kp)/2m* 
where m* — 1.13 m. Assuming long-lived quasiparticles, 
we approximate A(k,uj) by 6{uj — uj^) and perform the 
integration over the trap to obtain A{k,Lo) given by Q: 



A{k,uj) 



48fc2 



A(k,uj - iJ.{r)/fi) 



huj— air) 



(4) 



where /i(r) is the local chemical potential at position 
r. In order to calculate the integrated spectral function 
A{k, uj) of a Fermi liquid, we replace the spectral func- 
tion by 5(uj — LOk), and perform the integral in Q. kp{r) 
is calculated from the equation of state of the unitary 
gas determined in The temperature is chosen at the 
onset of superfluidity UbT/ijP = 0.32 0, |3i|. In order 
to make a direct comparison with the experimental data, 
we finally convolve our result with the experimental reso- 
lution in UJ 0, equal to 0.25Ep/h and results for various 
values of k are shown in Fig. 2] 

With no free parameter in the theory, FLT well re- 
produces the experimental spectra for A{k, lo) in the re- 
gion k < kp, with an excellent agreement in the re- 
gion 0.3feF < k < kp close to the most probable Fermi 
level in the trap (~ Q.lkp) where FLT is expected to be 
more accurate. Interestingly, we observe that the width 
of the peak at k/kp — 0.6 is well reproduced by our 
model meaning that the broadening of the line is not 
limited by the lifetime of the quasiparticles, but rather 
by trap inhomogeneity and measurement resolution. Sig- 
nificant deviations between experiment and FLT appear 
for k > lAkp, far from the most probable Fermi wave- 
vector. However in this region the energy spectrum signal 
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FIG. 4: Energy distribution data A{k,uj) from (red dots), compared with the prediction of Fermi liquid theory (blue lines), 
for k/kp = 0.3,0.6,0.9, 1.2, 1.5. 



is very broad and weak, corresponding to an incoherent 
background in the spectral function. Our Fermi liquid 
description thus accounts for the coherent part of the 
excitation spectrum from Q. 

In conclusion we have shown that the magnetic and 
thermal responses of the unitary Fermi gas support a de- 
scription of the normal phase in terms of Fermi liquid 
theory despite the fact that this system exhibits a high 
critical temperature for superfluidity. This behaviour is 
in contrast with underdoped cuprate high Tc materials 
displaying anomalous magnetic susceptibility or pseudo- 
gap physics in the normal phase. Recent quantum oscil- 
lation experiments on cuprates in high magnetic fields, 
aiming at studying the incipient normal state (somewhat 
analogously to the present work) do suggest long-lived 
quasiparticles [s^]. The drop of the susceptibility on the 
BEG side of the resonance for l/kpa ^ 0.5 indicates the 
appearance of a spin gap in this regime that deserves fur- 
ther investigations. Finally, the magnetic susceptibility 
could be a key observable for characterizing the onset of 
itinerant ferromagnetism in a repulsive Fermi gas [2^,[3l|. 
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AUXILIARY MATERIAL 

Measurement of the equation of state We recall 
here the procedure used to measure the equation of state 
P(/xi,/Z2,a), that was already employed in [l[ after the 
method proposed by [3| ■ We prepare a deeply degenerate 
mixture of the two lowest internal states of ^Li, held in a 
cylindrically symmetric hybrid optical/magnetic trap, of 
radial (axial) frequency iOr (<^z, respectively). The bias 
magnetic field Bq is chosen between 822 G and 981 G, 
allowing to tune the strength of interactions. The final 
atom number is 2 to 10 x 10'' atoms per spin state, and the 
gas temperature is 0.03(3)Tf, as measured from the fully- 
polarized wings of a trapped gas As shown in [ll,[3,|j], 
the local gas pressure along the z axis can directly be ob- 
tained from its in situ image. In the framework of local 
density approximation, this provides the grand-canonical 
equation of state P(/ii, /^2, a) at the local chemical poten- 
tials jiiz = — ^muj'^z'^, where is the global chemical 
potential for species i. The global chemical potential 
Hi for the majority species is directly obtained from the 
Thomas- Fermi radius Ri of the fully polarized phase, ac- 
cording to III = \mu)lR\. Similarly to [Ij], we obtain 
the global chemical potential by imposing that, at 
the outer radius i?2 of the minority species, the chemi- 
cal potential ratio uoZ/ui is given by the resolution of the 
impurity problem |8l4l3|. 

Fixed-Node Monte Carlo simulation The Haniilto- 
nian of the iVi + iV2 atoms of the two species is 

given by 

*2 / N^ N2 \ 

Ev? + Evn+En^-), (5) 

\i=l i' = l / 

where i, j, ... and i', j', ... label, respectively, majority and 
minority fermions. We model interspecies interatomic 
interactions using an attractive square well potential: 
V{r) ~ ~Vo if r < i?o and zero otherwise (Vq > 0). The 
short range i?o is fixed by the condition nR^ = 10~^, 
where n — ni -\- n2 is the total atom density. The depth 
Vq is instead chosen as to give the proper value of the 
scattering length a along the BEC-BCS crossover. We 
consider a system with fixed total number of particles 
[N — 66) in a fixed volume V — . 

Finite-size effects have been reduced using the tech- 
nique described in 

Gap and spin susceptibility Let us consider a system 
containing iVi spin up and N2 spin down particles. We 
define M ^ Ni- N2 and N = Ni + N2 the polarization 
and the total atom number and we note E{N, M) the 
energy of the system. If one assumes that the energy 
can be expanded in M then by symmetry the linear term 
vanishes and one gets E{N, M) = E{N, Q) + M'^ /2x + ■■■■ 
With this definition, x is then the spin susceptibility of 
the system. Indeed, adding a magnetic field h contributes 



to a —hM term to the energy and we immediately see 
that the energy minimum is shifted from M = to M = 
Xh. 

This argument is no longer true in the case of a gapped 
system. Indeed, polarizing a spin balanced system costs 
the binding energy of the broken pairs. This definition 
applies to any system composed of spin-singlet dimers, 
from a fermionic superfiuid composed of Gooper pairs, 
or a pure gas of uncondensed molecules, and leads to the 
following leading order expansion 

E{N, M) = E{N, 0) + |Af I A + ... 

To evaluate the spin susceptibility, we add as above a 
magnetic field h changing the energy into E — hM. We 
see that for /i 7^ 0, the potential is tilted but the energy 
minimum stays located at Af = (as long as \h\ < A 
corresponding to the Pauli limit pointed out by Glogston 
and Ghandrasekhar in the case of superconductors 0, 

El). 
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FIG. Al: Dependence of energy E with spin population im- 
balance M for gapless (top) and gapped (bottom) systems. 
Top. Full line: The dependence with spin imbalance is 
quadratic and the curvature is equal to the inverse of the 
spin susceptibility x- Dashed line: In the presence of a spin 
polarizing field h, the energy minimum is shifted to M = xh- 
Bottom: gapped system. Full line: Energy in the absence of 
external spin polarizing field. The slope is equal to the gap 
A. Dashed line: in the presence of a spin polarizing field, the 
energy profile is tilted but the minimum remains located at 
M = 0. 

Thermodynamic signature of the pseudogap The 

pseudogap phenomenon can be defined as a dip in the 
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density of state p{e) close to Ep reminiscent of the true 
cancellation of p inside the superfluid gap. We note A* 
the width of the dip, and for the sake of simplicity we 
assume that A* ^ UbTf- In a simple model where one 
assumes that the excitations of the system are described 
by the Fermi-Dirac distribution, one can show that the 
spin susceptibility of the system is equal to the density of 
state at the Fermi level. For small imbalances, the Fermi 
levels of the two spin states lie within the dip. The spin 
susceptibility is thus x — Po- When the imbalance is 
larger, the two Fermi surfaces are outside the dip (when 
Epi — Ef2 >• A*), and the pseudogap excitations do 
not contribute anymore to the thermodynamics of the 
system. In this case, the spin susceptibility is given by 
X = Po > p5 (see Fig. EH). 



Po 



e 



=1. 




FIG. A2: Top: sketch of the density of state in a pseudogap 
model. A footprint of the molecular state appears as a dip of 
width A* and depth po — Po in the density of state. Bottom: 
Polarization p as a function of the magnetic field b, with po = 
0.5po and A* = O.l/i. At low imbalance, the Fermi levels 
of the two spin species lie inside the dip. This results in 
a depletion of spin excitations and a reduction of the spin 
susceptibility. 

Polarizability of the unitary Fermi gas It is inter- 
esting to express our data for the unitary Fermi gas in 
the usual variables of condensed matter physics, namely 
the polarization p = (jii — n2)/{ni +712) as a function 
of (pi — p2)/2i?F, where Ep is the Fermi energy. These 
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FIG. A3: Polarization p of the unitary Fermi g func- 
tion of the magnetic field /ii — p2, normalized to the Fermi 
energy Ep ~ 7j.^/2m(37r^n)^/''. The red line corresponds to 
the relation p = |x(pi ~ P2)/2i5F, with x = 0.54. 

quantities are calculated from the thermodynamic func- 
tion at unitarity ho[b) = h{5 = 0, 6) according to 

_ /io(b) pi - P2 h 

lh^{h)~hh',{hy 2Ep ^ {hoib)^lbh'„ib)y/^- 

This requires to take the derivative of our experimen- 
tal data, which decreases the signal-to-noise ratio. We 
obtain the data plotted in Fig |A31 In the superfluid 
phase the polarization remains equal to and jump to 
p ~ 0.4 at the superfluid/normal transition (for pi — 112 — 
0.4 • 2Ep). The polarization then increases linearly with 
the magnetic field over a large polarization range, accord- 
ing to p= |x(pi - P2)/2Ef, with X = 0.54. 
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